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I.  INTRODUCTION 


The  purpose  of  this  report  is  to  illustrate  the  application  of  the 
K-pulse  concept  to  a  class  of  distributed-parameter  systems  which  can  be 
modelled  by  finite  lengths  of  non-uniform  transmission  lines.  The 
K-pulse  of  such  a  system  is  the  excitation  (input)  waveform  of  finite 
duration  which  yields  response  waveforms  of  finite  duration  at  all 
points  of  the  system. 

Numerical  techniques  using  finite  element  methods  ace  developed  to 
derive  accurate  approximations  of  the  K-pulse  and  response  waveforms  for 
uniform  and  non-uniform  transmission  lines.  Comparison  is  made  with 
exact  results,  where  these  can  be  obtained  using  other  methods,  to 
illustrate  the  accuracy  and  utility  of  the  method. 


II.  THE  R-MATRIX  OF  A  TWO-PORT  LINEAR  SYSTEM 

Let  a  two-port  be  represented  symbolically  as  in  Fig.  1,  denoting 
traveling  wave  amplitudes  (voltages)  at  ports  1  and  2  by  ai  (inward 
traveling  wave)  and  b-j  (outward  traveling  wave).  Then  the  scattering 
matri x  (S)  is  defined: 

b  =  (S)  a  ,  (1(a)) 
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Figure  1.  A  linear  two-port  system. 


When  cascading  two-ports,  a  more  useful  relation  is  given  hy 


(?) 


a  matrix  which  relates  the  right-  and  left-traveling  wave  amplitudes  at 
port  1  to  the  corresponding  right-  and  left-traveling  amplitudes  at  port 
2.  If  we  denote  this  matrix  by  (I12),  we  see  that  the  matrix  (Rin) 
relating  the  wave  amplitudes  at  the  first  port  to  those  at  the  nth  port 
is  just  the  product  of  n-1  matrices: 


(*ln)  =  (*12)  (*23)  •  •  •  (*n-l ,n)  . 


(3) 
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III.  FINITE  ELEMENT  ANALYSIS  OF  LINE  WITH  DISTRIBUTED  SHUNT  CONDUCTANCE 

It  is  customary  to  analyze  distributed  parameter  networks  by 
utilizing  lumped  constant  elements,  such  as  breaking  a  transmission  line 
into  lumped  element  T-  or  Pi-sections.  However,  we  shall  use  delay 
elements  corresponding  to  infinitesimal  lengths  of  transmission  line 
with  fixed  delay,  which  permits  the  final  R-matrix  of  the  system  to  be 
expressed  in  terms  of  polynomials  in  the  variable  z  =  exp(-2sx),  where  t 
is  the  element  delay. 

Referring  to  Fig.  2,  let  the  continuous  shunt  loading  be  modeled  by 
N  sections  of  loss-less  line,  each  of  length  aL  =  L/N,  with  a  lumped 
shunt  conductance  Gn  for  the  nth  section.  We  shall  adopt  the  symbol  dn 
to  indicate  the  amplitude  (voltage)  of  the  wave  traveling  to  the  right 
(dextra)  at  the  left  end  of  the  nt(l  section,  and  the  symbol  sn  to  denote 
the  amplitude  of  the  wave  traveling  to  the  left  (sinistra)  at  the  same 
ref  ence  plane.  For  a  typical  section: 
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and  Pn  =  Gn(Gn  +  2Y0)-1.  The  relation  between  input  and  output  planes 
of  Fig.  2  now  becomes 


(ID 


(R) 


a  z 


-N/2 


"Pq  ( 1 ~2po ) 


(R) 


(12) 


where 


N 

a  =  TT  ( l+wn) ,  and 
n=0 


l 


(R) 


( Rl )  (R2)  •  •  •  (^n) 


/  Rll(z) 

\  zR21(z) 


Rl2(z) 

zR22(z) 


(13) 


and  the  Rij(z)  are  polynomials  in  z  of  order  N-l. 
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pll(z)  =  Rn(z)  +  p0zR21(z) 


Pl2(z)  =  «i2(z)  +  PozR22(z) 

(15) 

p2l(z)  =  'PoRll(z)  +  (l-2p0)R2i(z) 

p22(z)  =  -P0r12(z)  +  (l-2p0)zR22(z) 

Hence  the  Pij(z)  are  all  polynomials  in  z  of  order  N.  Each  represents 
the  Laplace  transform  of  a  train  of  N  +  1  equally-spaced  pulses  with  a 
fixed  net  duration  T  =  2N  t  =  2L/C,  independent  of  N;  where  C  is  the 
wave  velocity  on  the  unloaded  line. 

Let  us  interpret  these  4  finite  duration  waveforms  for  the 
distributed  parameter  system,  in  the  limit  as  N  +  «. 

From  Eqs.  10  and  13,  it  is  easily  shown  that  if  s^+j  =  0,  or  the 
exiting  line  in  Fig.  2  is  terminated  in  the  characteristic  impedance  Z( 
of  the  unloaded  line,  Pn(z)  is  the  transform  of  that  special  input 
waveform  Pn(t)  of  finite  duration  applied  at  the  left  end  of  the 


loaded  line  which  produces  (a)  a  single  attenuated  (by  a  factor  1/a)  and 
delayed  (by  T/2  =  L/c)  impulse  at  the  exiting  terminal,  and  (b)  a  left 
reflected  wave  p2i(t)  with  finite  duration.  If  Pn(t)  is  applied  at  the 
right  end  of  the  line  of  Fig.  2  with  the  left  end  terminated  in  Z0,  then 
(a)  the  exiting  wave  at  the  left  is  the  same  attenuated  and  delayed 
impulse  as  above,  while  (b)  the  right  reflected  wave  produced  is 
-p12(t),  also  of  finite  duration.  We  thus  define: 


pll(t) 

.  t~' 

(Pll(z)} 

f  _1 

p2l(t) 

=  Jsj 

{P2l(z)> 

-p12(t) 

II 

I— * 

{p12U)> 

r1 

p2?(t) 

- 1 

(P22(z)> 

K-pulse  of  network 
Vi  pulse  of  network 
Tr  pulse  of  network 

/i 


(16) 


pulse  (K-pulse  under 
time  reversal ) 


The  property  suggested  for  p22(t)  means  that  the  time-reversed  from  or 
p22(T-t)  is  the  K-pulse  for  the  same  line  with  all  dissipative  elements 
replaced  by  negative  equivalents,  i.e.,  exchanging  -Gn  f0r  Gn.  In  a 
lossless  system,  we  shall  find  that  Pn(t)  =  P22(T-t). 

While  the  properties  above  hold  for  the  pij(t)  of  the  line-lumped 
approximate  of  any  order  N,  we  are  interested  in  the  limiting  form  of 
these  as  N  +  «,  and  how  accurately  this  limit  may  be  extrapolated  from  a 
finite  element  model  with  N  less  than  50.  These  questions  will  be 


addressed  in  Section  6,  where  computational  results  for  various  N  are 
compared  with  exact  waveforms,  derived  for  the  uniformly  loaded  line  in 
Appendix  II. 

Values  of  the  shunt  loads  Gn ( N )  must  be  determined  before 
calculations  for  the  finite  element  model  are  initiated.  These  values 
are  derived  in  Appendix  I  for  the  line  with  uniform  shunt  conductance  as 
well  as  the  line  with  a  linear  taper  of  shunt  conductance.  For 
simplicity,  we  have  to  this  point  assumed  entering  and  exiting  lines  in 
Fig.  2  have  the  same  characteristic  impedance  as  the  unloaded  line;  the 
case  where  these  differ  can  easily  be  accommodated  by  adding  at  most  two 
simple  R-matrices  for  junctions  between  dissimilar  lines  to  the  chain 
product  in  Eq.  12.  This  case  will  be  discussed  fully  in  Section  5. 


IV.  FINITE  ELEMENT  ANALYSIS  OF  LINE  WITH  NON-UNIFORM  CHARACTERISTIC 

IMPEDANCE 

A  line  with  continuously  varying  characteristic  impedance  is 
modeled  by  N  uniform  line  segments  in  cascade,  the  characteristic 
impedances  of  the  discrete  elements  matching  that  of  the  line  as  a 
staircase  function.  If  the  phase  velocity  is  non-uniform  as  well,  the  N 
sections  will  be  of  dissimilar  lengths,  but  constant  phase  delay.  We 
shall  consider  here  the  special  case  arising  when  a  non-uniform  line  is 
used  to  model  transmission  and  reflection  by  a  dielectric  slab  with  a 
varying  dielectric  constant.  Referring  to  Fig.  3,  the  non-uniform  line 
is  represented  by  N  uniform  sections  with  Zn,  3n  the  characterise 
impedance  and  phase  constant  of  the  n^  section.  The  lengths  ( aL )n  of 
the  sections  are  chosen  such  that  0n  (AL)n  *  (?aL  =  constant.  Thus,  the 
propagation  time  T/2  through  the  non-uniform  slab  is  composed  of  N  equal 
delays  t  for  each  section  where  T  *  2Nt. 

The  R  matrix  of  the  n™  element  of  Fig.  3  is  given  by 
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■  iv v,  'vvrrrr »■ 


where  rn  =  (Zn+1  -  Zn)/(Zn+1  +  Zn),  kn  =  (1  +  Zn/Zn+1)/2.  As  before 
(see  Eqs .  10-14), 


(20) 


where 


(R) 


(21) 


N 

and  a  =  II  kn 
n=0 


(R) 


(Ri)  (R2)  •  •  •  (Rn) 


/  Rn(z) 

\  zR2i (z ) 


Rl2(z)  \ 

zR22(z)  / 


(22) 


Since  the  R -j j ( z )  are  polynomials  in  z  of  order  N-l, 
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(24) 


Pll(z)  =  Rn(z)  +  r0zR2i(z) 

Pl2(z)  =  R12(z)  +  r0zR22(z) 

P21(z)  =  r0Rn(z)  +  zR21(z) 

P22(z)  =  TqRj2(z)  +  zR22(z) 

As  in  the  previous  section,  the  P-jj(z)  are  all  polynomials  in  z  of  order 
N;  the  interpretation  of  the  4  finite  duration  waveforms,  each  of  fixed 
net  duration  2Nt  =  T,  which  are  the  inverse  Laplace  transforms  of  the 
Pjj(z)  follows  as  in  Section  3,  preceding. 

The  values  of  rn(N)  and  kn(N)  must  be  determined  from  the 
particular  dielectric  constant  variation  of  interest  before  calculations 
of  the  Nth  order  approximants  are  initiated.  These  values  are  derived 
in  Appendix  II  for  the  line  representing  linear  variation  of  en  in  a 
dielectric  slab.  In  Section  6,  the  convergence  of  the  waveforms  as  N 
increases  will  be  examined. 


V.  TERMINATION  EFFECTS 


In  Section  3  and  4,  the  discrete  element  model  was  employed  to 
approximate  continuous  variation  of  conductivity  or  impedance  of  a 
non-uniform  line.  Discontinuous  changes  in  line  parameters  such  may 
arise  at  exiting  and  entering  ports  were  not  considered.  In  Section  3, 
it  was  assumed  that  the  line  with  distributed  conductance  was  inserted 
between  lines  of  the  same  characteristic  impedance  with  unloaded  line. 
The  resulting  R  matrix  of  the  system  has  unit  determinant.  If  entering 
and  exiting  lines  have  characteristic  impedance  Zs  and  Zd  respectively, 
the  R  matrix  for  the  system  is  given  by  (Rj). 


where 

ks  «  (l  +  z0/zs)/2,  rs  =  (z0  -  ZS)/(Z0  +  zs)  , 
kd  =  (1  +  Zd/Z0)/2  and 

rd  «  (zd  -  z0)/(zd  +  z0)  . 

When  Zs  =  Zd,  the  determinant  of  the  matrix  (Ry)  has  unit  determinant, 
property  which  holds  whenever  entering  and  exiting  lines  have  a  common 
characteristic  impedance.  In  Section  4,  the  characteristic  impedance 


of  entering  and  exiting  lines  were  not  generally  the  same,  so  that  the 

resultant  R  matrix  of  Eq.  2  does  not  have  unit  determinant,  but  a  value 

of  ZN+1  .  if  the  limiting  input  and  exiting  impedances  of  the  line  of 
~ ^o 

Fig.  3  are  Zo  and  Zn+1,  respectively,  then  we  can  obtain  the  Ry  for  new 

entering  and  exiting  impedances  Zs,  Zd,  respectively,  by  Eq.  24  where 

kS  =  (i  +  Vzs)/2.  rs  =  (Z0  -  ZS)/(Z0  +  Zs), 

kd  =  (!  +  Zd/ZN+1 ) /2 ,  and  =  (Z<j  -  Zm+^ ) / ( Z^  +  ZN+1). 

The  resultant  (Ry)  will  have  a  determinant  of  Zj/Zg.  If  Z^  =  Zs, 
entering  and  exiting  lines  have  common  characteristic  impedance,  then 
the  resulting  (Ry)  will  have  unit  determinant. 

With  reference  to  Fig.  4,  let  the  R-matrix  of  the  distributed 
parameter  network  referred  to  entering  and  exiting  line  characteristic 
impedance  Zs  and  Z^  be  derived  as  outlined  in  Eq.  24.  Then  the  K-pulse 
of  the  system  with  load  Z<j  and  generator  impedance  Zs  is  given  as 


K(t)  {pll(z)  +  rdPi2(z)  +  rsP2i(z)  +  rdrsP22(z)}  where  Td  = 

(zd  -Z L)/(Zd  +  Zl)  and  r$  =  (Z|_  -Zs)/(Z|_  ►  Zs).  Similar  expressions  can 
be  obtained  for  other  response  waveforms.  Numerical  results  of  K-pulse 
and  reflected  pulses  for  a  finite  transmission  lines  inserted  between 
lines  with  various  Zs  and  Zd  will  be  illustrated  in  Sec.  6. 


VI.  EXAMPLES 


To  illustrate  simple  examples  of  the  K-pulse  in  a  distributed 
parameter  system,  we  model  the  reflection  of  a  plane  wave  at  normal 
incidence  to  a  planar  dielectric  slab  by  a  network  in  which  a  length  L 
of  line  with  conductance  G  is  inserted  between  two  semi -inf i nite  lines 
with  characteristic  impedance  Zs  and  Zd.  The  K-pulse  and  reflected 
waveforms  are  then  calculated  using  a  N-section  line-lumped 
approximation  to  the  finite  transmission  line.  In  the  following, 
various  numerical  results  obtained  for  the  approximations  of  the  K-pulse 
and  response  waveforms  for  uniform  and  nonuniform  lines  are  presented. 
Comparison  is  made  with  exact  results,  where  these  can  be  obtained  using 
other  methods,  to  illustrate  the  accuracy  and  utility  of  the  methods. 

A.  UNIFORM  CONDUCTANCE  G 

Consider  first  a  lossless  line  with  length  L  shorted  at  one  end 
and  appended  to  a  semi-infinite  line  of  twice  the  intrinsic  impedance, 
corresponding  to  a  slab  with  dielectric  constant  er  =  4.  Thus,  G  =  0, 

7-s  =  2Zo  and  Zd  =  0,  where  Zo  is  the  intrinsic  impedance  of  the  lossless 
slab. 

In  Fig.  5,  the  K-pulse  consists  of  a  unit  impulse  followed  by  a 
second  impulse  of  1/3  with  a  delay  of  ?L/C,  corresponding  to  the  transit 
time  down  and  back  the  shorted  line.  The  reflected  waveform  is  the  time 
reversed  negative  of  the  K-pulse.  If  uniform  shunt  conductance  loading 


G  is  now  introduced  along  with  shorted  line,  such  that  the  total 
conductance  in  the  finite  line  equals  the  surge  admittance,  the  K-pulse 
and  reflected  pulse  are  shown  in  Fig.  6.  As  in  the  lossless  case,  a 
unit  impulse  begins  the  input.  Thereafter,  signals  returning  to  the 
input  junction  and  reflected  down  the  line  segment  are  cancelled  or 
"killed"  by  subsequent  input  contributions  until  the  final  signal 
reflected  from  the  shorted  end  is  returned.  Since  no  further  signal 
travels  down  the  line  segment,  the  shorted  section  is  at  rest  after 
At  =  2L/C.  The  K-pulse  or  "ki 1 1 -pulse,"  has  a  duration  equal  to  the 
round  trip  time. 

The  waveforms  of  Fig.  6  are  calculated  using  a  N-section 
line-lumped  approximation  to  the  uniform  line.  For  comparison,  the 
exact  results  are  also  presented  in  Fig.  6.  It  is  observed  tha-  me 
K-pulse  and  reflected  waveform  converge  rapidly  to  the  exact  res  T’s  *■ 
N  increases  from  N  =  10  to  N  =  40. 


Fig.  7  and  Fig.  8  present  the  K-pulse  and  reflected  pulse  for  a  uniform 
lossy  line  inserted  between  two  semi -infinite  lines  of  different 
characteristic  impedance,  corresponding  to  a  planar  slab  separating  two 
half  spaces  with  characteristic  impedances  Zs  and  Z<j  respectively.  Fig. 
7  shows  the  results  for  the  symmetrical  case  (Zs  =  Z<i  =  2Zo),  where  Z0 
is  the  intrinsic  impedance  of  the  unloaded  line.  It  is  again  observed 
that  as  N  increases  from  1(1  to  40,  the  N-section  line-lumped 
approximation  rapidly  converges  to  the  exact  results.  Note  that  for  the 
symmetrical  case,  the  left-reflected  and  right-reflected  pulse  are 
identical.  For  the  asymmetrical  case  shown  in  Fig.  8  (Zs  =  2Zo,  Zd  = 

Z0),  though,  the  left  and  right-reflected  waveforms  are  not  the  same. 

In  both  cases,  the  K-pulse  and  response  waveforms  have  a  duration  time 
equal  to  the  round-trip  transit  time. 


EXACT  SOLUTION 


Figure  7a.  K-pulse  and  reflected  waveforms  for  uniform  lossy  slab 
( Zs  =  Zd  =  2.0,  k  =  1) 


Figure  8c.  K-pulse  and  reflected  waveforms  for  uniform  lossy  slab 
Us  =  2.0,  Zd  =  0.5,  k  =  1) 


B.  LINEARLY  VARYING  CONDUCTANCE  G 

Consider  a  finite  line  with  linearly  tapered  shunt  conductance, 
with  the  total  shunt  conductance  equal  to  the  intrinsic  impedance  of  the 
unloaded  line.  Again,  the  continuous  shunt  loading  is  modeled  by  N 
sections  of  lossless  line,  each  of  length  M  =  L/N,  with  a  lumped  shunt 
conductance  Gn  for  the  nth  section.  Values  of  the  shunt  loads  Gn  are 
determined  in  Appendix  II.  Figs.  9-11  present  the  K-pulse  and  response 
waveforms  for  a  linearly  tapered  line  with  various  terminations.  It  is 
observed  that  convergence  to  the  K-pulse  and  response  waveforms  are 
rapid  and  simple,  even  for  lines  with  continuously  varying  parameters. 
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.  K-pulse  and  reflected  waveforms  for  shorted  tapered  line 

(Zs  =  2.0,  Zd  =  0.0,  k  =  1) 


C.  LINEAR  VARIATION  OF  CHARACTERISTIC  IMPEDANCE 


i 
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Numerical  results  of  K-pulse  and  reflected  waveforms  for  a  line 
with  continuously  varying  characteristic  impedance  are  presented  in  this 
section.  Waveforms  are  obtained  using  a  model  with  N  uniform  line 
segments  in  cascade,  the  characteristic  impedances  of  the  discrete 
elements  matching  that  of  the  line  as  a  staircase  function.  If  the 
phase  velocity  is  non-uniform  as  well,  the  N  sections  will  be  of 
dissimilar  lengths,  but  constant  phase  delay.  We  shall  consider  here 
the  special  case  arising  when  a  non-uniform  line  is  used  to  model 
transmission  and  reflection  by  a  planar  dielectric  slab  with  a  linearly 
varying  dielectric  constant.  The  profile  of  the  dielectric  constant  er 
is  shown  in  Fig.  12.  The  two  half  spaces  separated  by  the  slab  have  a 
dielectric  constant  es  and  respectively. 


Fig.  13  shows  the  K-pulse  and  reflected  pulse  for  a  planar  slab 
with  a  dielectric  constant  er,  linearly  tapered  from  es  =  eQ  to 
ed  =  2 e0.  The  waveforms  shown  in  Fig.  13  are  calculated  using  a 
20-section  line-lumped  approximation  to  the  planar  slab.  Calculations 
using  N  =  40  shows  that  convergence  to  the  K-pulse  and  reflected 
waveform  is  rapid  and  simple. 


VII.  CONCLUSIONS 


In  this  report  we  illustrate  the  application  of  the  K-pulse  concept 
to  a  class  of  distributed-parameter  systems  which  can  be  modelled  by 
finite  lengths  of  non-uniform  transmission  lines.  The  K-pulse  of  such 
a  system  is  the  excitation  (input)  waveform  of  finite  duration  which 
yields  response  waveforms  of  finite  duration  at  all  points  of  the 
system.  Numerical  techniques  using  finite  element  methods  are  developed 
to  derive  accurate  approximations  of  the  K-pulse  and  response  waveforms 
for  uniform  and  non-uniform  transmission  lines.  Comparison  is  made  with 
exact  results  derived  for  the  uniformly  loaded  line,  to  illustrate  the 
accuracy  and  utility  of  the  method. 

The  next  logical  step  in  this  analysis  is  to  address  the  inverse 
problem,  i.e.,  given  the  K-pulse  and  response  waveforms,  what  are  the 
electrical  parameters  of  the  line.  Kennaugh  had  claimed  in  an  earlier 
report  [4 1  that  synthesizing  the  parameters  of  the  non-uniform  line  from 
measured  K-pulse  and  response  waveforms  was  as  equally  tractable  as  the 
direct  problem.  The  key  role  of  the  K-pulse  in  factoring  the  system 
before  attempting  synthesis  has  been  clearly  established  in  this 
approach,  which  differs  from  the  one-dimensional  inversion  techniques. 

It  is  intended  to  investigate  this  problem  as  time  and  funds  permit. 


APPENDIX  I 


i 

K-pulse  for  Uniform  Lossy  Line  with  Short-Circuit  Termination. 

Reflection  by  a  lossy  distributed-parameter  network  furnishes  a 
**  useful  example  for  application  of  K-pulse  concepts.  As  shown  in  Fig 

A-l,  a  length  1  of  lossy  line  with  characteristic  impedance  Z-|  is 
shorted  at  the  far  end.  The  reflection  coefficient  (voltage)  at  the 
input  terminal  when  connected  to  a  uniform  line  of  characteristic 
impedance  Z0  is  of  interest. 


SHORT 

CIRCUIT 
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Figure  A-l.  Lossy  line  configuration 
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We  assume  that  line  loss  is  modeled  by  a  uniformly  distributed 
shunt  conductance,  such  that  the  total  shunt  conductance  over  the  length 
1  is  a  specified  fraction  of  the  surge  admittance  of  the  line  without 
loss,  i.e. 

total  Gshunt  =  k 

n 

The  ratio  of  the  surge  admittance  of  the  entire  line  (to  the  left 
of  reference  plane!  to  that  of  the  line  is  also  a  specified  parameter: 


The  expression  for  the  voltage  reflection  coefficient  can  be  derived 
as: 


/  K  +  St 

/K  +  St 

rv  =  u  tanh 

(St 

VS)- 

Y  S 

/  K  +  St 

/  K  +  St 

u  tanh 

(St  \ 

/  St  )  + 

V  St 

l 

where  S  =  jurr,  t  =  TTjf  and  Ci  is  the  velocity  of  a  wave  traveling  the 
finite  line  section  in  the  absence  of  loss. 


Recognizing  that 


ez  -  e-z 

tanh  (*)  =  ^  T~z 


the  voltage  reflection  coefficient  rv  can  be  rewritten  as 


- -  ./St+K  - -  ^  /St+K 
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Multiply  numerator  and  denominator  by 

-S  T 

e  ,  one  obtains 
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with  a  change  of  variable,  Y  =  St  +  a,  where  a  =  K/2,  rv  can  be  written 


It  can  be  seen  that  the  expressions  in  the  numerator  and  denominator  of 
the  expression  given  for  rv  are  entire  functions  of  S.  The  inverse 
lapalace  transform  of  the  denominator  gives  the  K-pulse  K(t)  while  the 
reflected  pulse  r(t)  is  the  inverse  transform  of  the  numerator,  i.e. 


t  - 

eatK(t) ^K(s)  *  e® 


Ry  using  the  tables  for  inverse  transforms: 
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obtain 


-  \K\J\2xJ  \2xJ  J  |_u(t)  -  u^~2x 


If  we  set  2t  =  1,  the  total  duration  of  R ( t ) 
the  normalized  K-pulse  is  given  by 


K (t )  =  S(t)  -  e“K 


5(t-l) 


-  2J0  (Kh/tT2  )  CU(t)  -  U(t-l)]  , 


APPENDIX  II 


ELEMENT  VALUES  FOR  DISTRIBUTED  SHUNT  CONDUCTANCE  LINES 


The  relation  between  the  R  matrix  and  S  matrix  has  been  indicated 
in  Sec.  2.  The  elements  of  the  S  matrix  for  the  finite  element  model  of 
the  distributed  line  is  given  by 

Yc  -  Yin  . 

S11  =  Yc  +  Yin 

sl?  =  $21  =  1+S11  » 


where  Yc  is  the  characteristic  admittance  of  the  line.  For  a  typical 
section  of  the  line  with  a  shunt  conductance  load  Gn>  Yin  =  Yc  +  Gn, 
thus 


Sn  =  -$n  =  -Gn/Yc 

7Tc~+l^  2  +  Gn/Yc 


In  the  case  of  a  line  with  distributed  shunt  conductance,  we  assume  that 
the  total  shunt  conductance  Ggh  equals  to  a  specified  fraction  of  the 
characteristic  admittance,  i.e.,  Gsh  =  kYc.  Then,  the  element  value  Gn 
for  a  N  section  model  is  simply: 

^n  =  kyc  for  uniform  lossy  line, 

N 

and 

p 

’n  =  ?kN  Yc  for  linearly  tapered  line. 

nTn+TT 

With  the  Gn  given  as  above,  it  is  straightforward  to  calculate  the 
R-matrix  elements  described  in  Sec.  3. 


APPENDIX  III 


ELEMENT  VALUES  FOR  VARIABLE  CHARACTERISTIC  IMPEDANCE  LINES 

The  approximate  K-pulse  and  response  waveforms  for  a  line  with 
continuously  varying  characteristic  impedance  are  derived  by  using  a 
model  with  N  uniform  line  segment  in  cascade,  the  characteristic 
impedances  of  the  discrete  elements  matching  that  of  the  line  as  a 
staircase  function.  For  the  special  case  arising  when  a  non-uniform 
line  is  used  to  model  transmission  and  reflection  by  a  planar  dielectric 
slab  with  a  linearly  varying  dielectric  constant  er,  the  dielectric 
constant  en  for  the  nth  section  of  the  N-section  model  is  simply 

en  =  (£2  -  £1)  (2n-l)  +  ei 


where  is  the  dielectric  constant  corresponding  to  the  entering  line, 
and  ^2  is  that  corresponding  to  the  exiting  line.  With  en  given  as 
above,  the  element  for  the  R  matrix  can  be  readily  determined  as 
described  in  Sec.  4. 


L 


i 


REFERENCES 

I 

(1)  "Principles  of  Microwave  Circuits",  Chap.  V,  Sec.  5.1ft,  Radiation 
Laboratory  Series,  8,  Montgomery,  Dicke  and  Purcell,  Ed., 
McGraw-Hill  Book  Company,  1948. 

M  (2)  "Acoustic,  Electromagnetic  and  Elastic  Wave  Scattering  -  Focus  on 

the  T-matrix  Approach",  Varadan  and  Varadan,  Ed.,  Pergamon  Press, 
1980. 

(3)  "Basic  Theory  of  Waveguide  Junctions  and  Introductory  Microwave 
Network  Analysis",  Kerns  and  Beatty,  Ed.,  Pergamon  Press. 

(4)  "Fifth  Annual  Report  710816-12",  December  1982,  The  Ohio  State 
University  ElectroScience  Laboratory,  Department  of  Electrical 
Engineering;  prepared  under  Contract  No.  N00014-78-C-0049  for  the 
Department  of  the  Navy,  Office  of  Naval  Research,  Arlington, 

t  Virginia  22217. 


II 


I 


END 


s 


FILMED 


r 


6-85 


DT1C 


’  V' 


